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An effective model with constituent quarks as fundamental degrees of freedom is used to predict 
the relative strangeness production pattern in both high energy elementary and heavy ion collisions. 
The basic picture is that of the statistical hadronization model, with hadronizing color-singlet clus- 
ters assumed to be at full chemical equilibrium at constituent quark level. Thus, by assuming that 
at least the ratio between strange and non-strange constituent quarks survives in the final hadrons, 
the apparent undersaturation of strange particle phase space observed in the data can be accounted 
for. In this framework, the enhancement of relative strangeness production in heavy ion collisions 
in comparison with elementary collisions is mainly owing to the excess of initial non-strange matter 
over antimatter and the so-called canonical suppression, namely the constraint of exact color and 
flavor conservation over small volumes. 



I. INTRODUCTION 



The recent observation that hadron multiplicities in e + e~ and hadronic high energy collisions agree very well with a 
statistical- thermodynamical ansatz |iJ.|2{ has revived the interest on the statistical models in high energy collisions, 
an idea dating back to the 50's 0,0. As the apparent chemical equilibrium of all hadron species in these elementary 
collisions (EC) cannot be driven by inelastic collisions amongst hadrons after their formation, this finding has led 
to the idea of a pure statistical filling of multi-hadronic phase space of homogeneous hadronizing regions {clusters 
or fireballs) as an intrinsic feature of hadronization process, occurring at a critical value of energy density HHH- 
Otherwise stated, hadrons are born in an equilibrium state, as envisaged by Hagedorn 7]. It must be emphasized 
that, within this framework, temperature and other thermodynamical quantities have an essential statistical meaning 
which does not imply the occurrence of a thermalization process via inelastic hadronic collisions. 

The same model, in different versions, has been successfully applied to a large set of heavy ion collisions (HICjdata 
H, 0, ^3 and this has been interpreted as a clue of the minor effect of post-hadronization inelastic rescattering H , an 
indication which is also supported by kinetic models calculations [Tl| . 

One of the nice features of the statistical hadronization model (SHM) is the very low number of free parameters 
required to reproduce a large number of hadronic multiplicities. Provided that the masses and charges of the assumed 
hadronizing clusters fluctuate according to a particular shape function pl.ll2j|. all Lorentz-invariant quantities (like e.g. 
hadron average multiplicities) depend on two parameters: the sum V of the volumes of the clusters and the temperature 
T . However, in order to get a satisfactory agreement with the data, the model has to be supplemented, both in EC 
and HIC, with one more phenomenological parameter, 75, suppressing the production of particles containing n strange 
valence quarks by 7^ with respect to the expected production in a fully equilibrated hadron gas . 

The behavior of 7s as a function of collision and center-of-mass energy was found to be rather odd as it turns out 
to be significantly higher in c + c~ collisions than in hadronic collisions over a large energy range 0,0 whereas it has a 
fairly stable value in heavy ion collisions 0, which is approximately the same as in e + e _ collisions. A clearer insight 
in the mechanism of strangeness production can be achieved by estimating the ratio between newly produced valence 
strange and up, down quark pairs, the so-called Wroblewski factor A5: 

2(ss) 

AS - <uu) + (dd) (1) 

which shows a striking regularity, being in fact fairly constant in all kinds of elementary collisions |lf)j | over two orders 
of magnitude of center-of-mass energy and as twice as large in high energy heavy ion collisions (see Fig.[^). It should 
be stressed that this ratio is calculated by counting the primarily produced quark pairs, i.e. those belonging to directly 
emitted hadrons. As primary hadrons are not measurable, they must be calculated by using a model and As turns 
out to be a model dependent quantity especially with regard to the number of u and d quarks which significantly 
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increases during the post-hadronization hadronic decay chain. However, if the number of measured species is large 
and the model accurately reproduces them (which is the case for the SHM) , the thereby estimated As is expected to be 
reasonably close to the actual value. The need of an extra suppression of strangeness with respect to the full statistical 
equilibrium for the hadronic system urges the search for an explanation based on a more microscopic approach. This 
is indeed the main subject of the present work, where we try to calculate As within one statistical model scheme 
in both EC and HIC, employing a constituent quark model as the basic underlying structure, with the purpose of 
justifying the lack of strangeness chemical equilibrium at hadron level observed in the data. In fact, constituent quark 
models have already been used to calculate strangeness production in HIC on the basis of transport equation |l3j . 
Amongst other microscopic approaches to this problem and, more generally, to account for hadronization equilibrium 
features, a recent study has been performed in ref. focussed on the role of massive Polyakov loops. 

The paper is organized as follows: in Sect. [H] we introduce the physical picture and in Sect. Hill the full model is 
described. In Sect. lIVl the Nambu-Jona-Lasinio model with exact conservation of quantum charges is used to calculate 
As and compare it with the corresponding value obtained from fits of the statistical hadronization model to the data 
in both HIC and EC. Also, the stability of the results is addressed. Sect.^is devoted to conclusions. 



II. THE PHYSICAL PICTURE 



The physical picture keeps the same scheme of the SHM: formation of a set of clusters endowed with charge, 
momentum, mass and volume, in local statistical equilibrium. However, statistical equilibrium is now assumed to 
apply to a system of constituent quarks whilst hadrons are assumed to be produced via their coalescence, still in a 
purely statistical fashion (statistical coalescence) so as to give rise to a partially chemically equilibrated hadron gas. 
Furthermore, each cluster is required to be a color singlet. This requirement has of course no effect at hadron level, 
but it is crucial in a quark model. 

Two more crucial assumptions are introduced. First, the thermodynamical parameters fitted in the SHM, in high 
energy EC and HIC, are interpreted as the critical values for color deconfinement and, consequently, for (approximate) 
chiral symmetry restoration [l5lll6j| : secondly, it is assumed that the produced s and light quarks, or at least the ratio 
s/u, i.e. As, survive in the final hadrons after hadronization has taken place. 

The first assumption is supported by the constancy of the fitted temperatures in the SHM for various processes in 
EC and by its agreement with other estimates of the critical tem per ature; in view of this fact, the use of effective 
models such as the Nambu-Jona-Lasinio model [181 . 1191 and others poll2lll22ll2^| embodying chiral symmetry breaking 
and restoration, looks well suited. It should be pointed out that the values of the thermodynamical parameters (i.e. 
temperature and baryon-chemical potentials) extracted within the SHM in HIC (see Table lie in a range where all 
effective models and recent lattice simulations |2^| predict smooth cross-over transitions, far from a possible critical 
ending point. For the same reason, extensions of the NJL model |26| at high baryon chemical potential and low 
temperature, where color superconductivity takes placej^j, can be disregarded. 

The second assumption allows us to calculate As within a quark model and compare it with that obtained from 
valence quark counting in the produced hadrons. Also, it implies that the lack of full chemical equilibrium at hadron 
level is indeed the consequence of full chemical equilibrium at constituent quark level. 

In the grand-canonical framework, which applies to a large system, As would depend only on intensive quantities 
such as chemical potentials and temperature and its constancy would be tightly related to theirs. However, if the 
system is not large, the canonical (and possibly micro-canonical) ensemble, in which the exact conservation of cluster 
quantum numbers in the hadronic or quark system is enforced, must be used and, as a consequence, As gets a 
dependence also on the volume (the so-called canonical suppression). While the volume within which quantum 
charges (baryon number, electric charge and strangeness) are fixed can be taken as the sum of the proper volumes of 
all clusters, though under appropriate non-trivial assumptions [2lll2|. the volume over which color must be neutralized 
is, as has been mentioned, that of a single hadronizing cluster. This very fact introduces a further parameter in the 
model, i.e. the average single cluster volume V c which can indeed be much smaller than the volume over which flavor 
is conserved. As quarks do carry color charge, As might be significantly affected by variations of V c . 

Therefore, the underlying idea is that the presence of a characteristic constant value of V c in EC, related to the 
typical distance over which color is neutralized, may sizeably reduce the value of As (together with flavor constraint) 
whereas color deconfinement in HIC may lead to an enhancement of s/u ratio. We thus argue that the main difference 
between hadronization in EC and HIC is to be found in the typical size of the pre-hadronization color-neutral region: 
it has to be something of the order of a hadron radius in EC (a sort of mini-QGP droplet) and much larger in 
HIC if a macroscopic (i.e. extending over many hadron volumes) QGP has been formed. Of course, it must be 
emphasized that the process stage at which the statistical equilibrium is achieved is expected to be deeply different in 
EC as compared with HIC: an early thcrmalization of partons is indeed envisaged in HIC, which is maintained until 
hadronization, whereas a late, pre-hadronization local equilibrium scenario, possibly driven by the strongly coupled 
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non-linear evolution of QCD fields in the late soft regime is envisaged for EC. 

Hence, we argue that relative strange quark production might probe these two scenarios because of a possible 
stronger color and flavor-canonical suppression of s quarks with respect to u and d quarks in EC which is absent 
in HIC due to much larger single-cluster volume and overall volume. The ultimate reason of the stronger canonical 
suppression of s quarks in comparison with u, d, resides in their different constituent mass values. The hereby proposed 
mechanism of strangeness enhancement in HIC is rather different from that suggested by Mullcr and Rafelski j2^| 
insofar as no difference in the time scale is invoked but an equilibrium situation is assumed in both EC and HIC, 
though over differently sized regions. 



III. THE MODEL 



According to the basic idea of the SHM and what has been said in the previous section, the formation of a set of pre- 
hadronic, color-singlet clusters at statistical equilibrium, with definite values of flavor quantum numbers is envisaged 
as the result of a complex dynamical evolution. In an equilibrium scheme, all physical observables relevant to a 
particular cluster can be calculated by means of suitable operations on its associated partition function. Provided 
that cluster masses are large enough so that micro-canonical effects can be neglected, the requirement of definite 
quantum numbers implies that the i th cluster's partition function has to be the canonical one rather than the more 
familiar grand-canonical: 

Zi= £ e- B/T (2) 

states 

where T is the cluster's temperature, E its energy and the sum is meant to run over all allowed states, namely the 
states with flavor and color quantum numbers matching those of the cluster. The calculation of Zi can be done by 
using the symmetry group G associated to the involved quantum numbers. If the continuous group G is an exact 
internal symmetry for the Hamiltonian H and the cluster state belongs to its irreducible representation v, then |2flj |: 

Zi= J dfi(ji, ...,7r)Xw(7i, -,7r)Zi(7i, ...,7 r ) = / d/i(7i,...,7 r )x^(7i,->7r)Tr [e J (3) 

where dfifyi, •■•,7 r ) is the invariant normalized measure of G, Xv is t Qe character of the representation v, r is the rank 
of G, Qi are the mutually commutating generators of the Cartan subalgebra and 7/ are the group parameters. The 
actual group of interest is SU(3) xU(l) 3 , where each U(l) corresponds to conserved net flavors U, D and S and SU(3) 
to the color. For instance, the canonical partition function of the i th cluster of a gas of free quarks (no antiquarks 
included) with flavor j involves the following Zi in Eq. © : 



Zi(Vi, #2, <Pj) = exp 

where: 



£^»,oK*)e'"*' 7 £/<f>' T ) (4) 

\7l— 1 



Xi,oO?i, 2 ) = + + e 1 ^ + °*) (5) 

is the character of the fundamental representation of SU(3), with 9\ 2 € [— 7r, it]; 4>j is the parameter of U(l)j 



V is the volume, g is the spin degeneracy factor, Cj — yp 2 + m^ and the sum over discrete states k has been 
approximated with its continuum limit. When including antiquarks and considering three flavors (i.e. the full group 
SU(3)xU u (l)xUd(l)xU s (l)), Z l can be written as the product of three functions like those in Eq. 

Zi (0i,e 2 , 4> u ,4> d , <j,„) = JJ Za {e 1 ,e 2 , ( f )j ) {Ox, e 2 , <t> 3 ) (6) 

j—u 7 d : s 

The irreducible representation v can be labelled by three integer numbers, the net flavor Ui,Di,Si, and, as far as 
the SU(3) color group is concerned, v is simply the singlet representation. Therefore, by introducing the vectors 
Q.; = (Ui, Di, Si) and 4> i — ((pur, 4>iD, 4>is), the character can be written as \v = ex P[ — l Qi ' 4>i\ an d: 



Z%j — 
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where: 



dn{p u 2 ) = 



3! 



d9i 
27 



^2 

2tt 



3 



sm' 



2 J 



(8) 



with the constraint Yl^—i Ok — mod 2-7T, is the SU(3) invariant integration. The partition function hitherto considered 
is relevant to one cluster. Yet, there are several clusters in a single collision event and, as long as global observables 
are concerned, theoretical predictions require summation over all clusters. Moreover, clusters may well be produced 
with different configurations of flavor numbers (Ui, Di, Si) in different events (provided that their sum fulfills the 
conservation law, i.e. it must be equal to the flavor numbers of the colliding particles), hence an integration over 
all possible configurations weighted by its probability must be carried out in order to calculate averages of physical 
quantities. However, the probability distribution of clusters flavor configuration is an unknown function, which cannot 
be predicted within the statistical model and which is most likely governed by the preceding dynamical process. We 
will therefore assume a configuration probability distribution w which is the "maximum disorder" one and which allows 
a remarkable simplification of the expressions of averages quark multiplicities as well as any other Lorentz-invariant 
observables. This function w, for an event with N clusters can be written as: 



rii=i ^(QOSq.EiCK 



«»(Qi 



Q 



N) 



E 



Qi,...,Qj» 



(9) 



where Zi as in Eq. J7J), the S ensures that global conservation of initial quantum numbers Q is fulfilled. If one has to 
calculate the average multiplicity of a given quark species j in a N cluster event, it is advantageous to introduce a 
fictitious fugacity Xj and taking the derivative of log Zi with respect to it, so that the overall multiplicity reads: 



d 



w(Q 1 ,...,Q N )J2—logZ i {Q i 

i—l 



Qi,-,Q 



V to(Q x ,...,Qjv) ^-TTlog^(Qi) 
o o 



N 



i=l 



A,-=l 



A,=l 



(10) 



and, by using Eq. ©: 



dX 



log 



N 

E 

Qi,...,Qjv »=1 



A 4 =l 



(11) 



Were not for the color-singlet constraint stuck to each Zi, the sum over all configurations on the right hand side 
of Eq. Qllj l would be, for Xj = 1, the canonical partition function of a single cluster having as volume the sum of 
volumes of the individual clusters and with quantum number vector Q ,2] , provided that all clusters have the same 
temperature. Thereby, also the explicit dependence on N would vanish. However, this is no longer true once the 
color singlet constraint is introduced. For a similar strong equivalence with one cluster to apply, single clusters should 
be in suitable color quantum states superpositions, which is definitely against the common belief of prc-confincmcnt 
and local color neutralization. Nevertheless, it can be proved that a considerable simplification in Eq. 1)11(1 occurs if 
clusters in one event are assumed to have the same volume, i.e. Vi = V c . In that case, it can be proved (see Appendix 
A) that the right hand side of Eq. (|ll|l . which can be defined as the global partition function, reads: 



Z = 



ny if 



d/*(0i,0a) Z i {9 1 ,e 2 ,ct>) 



v/v a 



(12) 



where V = J^- Vi = NV C (for the free quarks gas Zi is given by Eq. © with the single flavor term Z i j{9 1 ,02,4>j) 
calculated from Eq. (0J with (J5J). The introduction of the local color-singlet constraint gets the value of the global 
partition function severely reduced with respect to an only globally color-singlet constrained case where clusters are 
allowed to be in a quantum superposition of colored states. This can be quite easily understood for any extra constraint 
essentially involves a decrease of the number of available states in the system. If V c = V and V is very large, i.e. 
in the thermodynamical limit, it can be shown that the usual grand-canonical expressions of the partition function, 
average multiplicities etc. are recovered. 
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Besides the basic assumption of local statistical equilibrium (for a single cluster), the expression (|12fl ultimately 
relies on a particular choice of flavor distribution (i.e. jS}) among clusters, which is indeed a non-trivial assumption. 
This particular distribution allows the reduction to one equivalent global cluster just because it is the "maximum 
disorder" distribution, namely the probability of obtaining a configuration (Qi, • • • , Qjv) by randomly splitting a 
statistically equilibrated global cluster into N sub-cluster each with volume V c and in a color singlet state ; in fact, 
this can be proved by showing that the w's in Eq. JHJ) minimize the free energy of the system (see Appendix B). 

It should pointed out that hadronizing clusters might be too small for a local (i.e. in each of them) temperature to be 
defined, as we have tacitly assumed. Should this be the case, the appropriate treatment would be micro-canonical and 
not canonical, with clusters described by mass and volume instead of temperature and volume. Notwithstanding, as 
far as the calculation of Lorentz-invariants is concerned, a global temperature may still be recovered by resorting to a 
similar equivalent global cluster reduction procedure at micro-canonical level, which is described in detail in ref. |l2j 2 , 
with a micro-canonical equivalent of the distribution @ . The proof is quite lengthy and the mathematical framework 
more involved, but it is no longer necessary to assume that each cluster ought to have the same temperature. In fact, 
only the equivalent global cluster or, equivalently, the global canonical partition function lJT2)l. must be large enough 
to allow a canonical treatment of the system as a whole. 

We now have to face the problem of a suitable choice of a microscopic Hamiltonian to be used in Eq. © and, 
thereby, make a theoretical calculation of As. As full QCD is out of question, we resort to QCD-inspired low energy 
models. The simplest model to start with is the Nambu-Jona-Lasinio (NJL) model, which has quarks as fundamental 
degrees of freedom and embodies essential features of dynamical chiral symmetry breaking (xSB) in the limit of 
vanishing current quark masses. Also, we confine ourselves to mean- field approximation, as presented in ref. jl9j . 
using the effective Hamiltonian: 

Hmfa = V [g s (a 2 + /3 2 + 7 2 ) + 4g D a/3 7 ] + J d 3 x q (-17 • V + M) q (13) 

where V is the spatial volume, g s and g D are the four-fermion and the six-fermion J7(l) J 4-breaking couplings respec- 
tively [HOH, a = (uu},(3 = (dd) , 7 = (ss) are the quark condensates and: 



M u = m u 


— 2g s a - 


- ^9 D h 


M d = m d 


-2g s P- 


2g D ui 


M s = m s 




2g D a(3 



(14) 

are the constituent quark masses which, owing to the contact four-fermion interaction, are linearly related to the 
current quark masses m u , m c i, m s . As far as the grand-canonical calculation is concerned, the variational parameters 
a, (3, 7 are determined by minimizing the effective potential ftQ] obtained within standard methods |3lll32l|. Of course, 
the same result can be readily obtained by calculating: 



Z = Tr e 



(15) 



with H = H MFA . In Eq. (|15|) . Nj are the conserved flavor operators and r| m i n is the effective action at the physical 
point. For constant fields r| min — f3VV\ m m = —logZ, where V is the effective potential. 

The quark over antiquark excess is obtained by taking the derivative of V with respect to the chemical potential for 
a given flavor: 

(nj - nj) = -— (16) 

whereas the number of particles, in the grand-canonical ensemble, reads: 

«-y f\ P , (1 „ 



(i 



/( V /P 2 + M|- M j) +1 



2 Therein, the proof has been given for a hadron gas but it can be easily extended to the present case. 
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where A is an UV cutoff and N c the number of colors. 

One of the remarkable features of the Hamiltonian H MFA in Eq. (|13[1 is that the previously obtained expressions in 
Eqs. Q, l|12|) for the free quark gas partition functions with exact flavor and color conservation, can be taken up by 
simply replacing the current masses with the constituent masses wij — » Mj and setting the cutoff A as upper bound 
for the momentum integration in Eq. l|4")l. 

The NJL model is the simplest choice in order to study strange quark production, yet it is not a compelling 
one. Actually, its validity is limited within a temperature range well below the UV cutoff A j2(|. However, we are 
reasonably confident that the main results found within this model would not essentially change when employing 
other effective models whose validity, in principle, e xtend to arbitrarily high temperatures. As an example, it is worth 
mentioning the model developed in refs. |2Ct I2H |22| . named in ref. as ladder-QCD. In this model the exchange of a 
gauge particle between quarks is considered, implying a momentum dependence of the self-energy and, consequently, 
an irreducibility to a simple model with free constituent quarks as in NJL. However, ladder-QCD is considerably 
different from NJL only in the UV regime, where a 1/p 2 tail is added to the self-energy, which is otherwise constant in 
the small momentum region |2(ll33| . The effect of the different behaviour in the UV regime on the integrated number 
of particles is expected to be rather small [20| . 



IV. ANALYSIS AND RESULTS 



The goal of the numerical analysis is the calculation of As in both HIC and EC, by performing the ratio between 
newly produced strange quark and u, d quark pairs in the NJL model. The main ingredient needed to calculate As are 
the constituent quark masses which, for a given temperature and baryon-chemical potential, have to be determined by 
the minimization of the free energy F and thus depend on the parameters contained in the effective Hamiltonian (TSJl, 
i.e. gs, A, go and the current quark masses. Their values have been fitted at zero temperature and baryon density 
in ref. |19| to many static and dynamical meson properties: 

m = m " ~^ md = 5.5 MeV; m s = 135.7 MeV 

A = 631.4 MeV; g s A 2 = 3.67; g D A 5 = -9.29 (18) 
Particularly for go , in ref. jl9j a temperature dependence is introduced through the phenomenological law: 



^(T) = 5D (T = 0)e-™ (19) 

where To is a further free parameter and gn(T = 0) is quoted in Eq. (JTSJ. Different values of To give rise to different 
behaviors of the constituent u, d mass (see Fig. plater on) as a function of the temperature as well as different positions 
of the cross-over curve in the (/jLb,T) plane. According to the statement in Sect.[Hl we take the thermodynamical 
parameters T and /is in Table [I] as an input, and therefore To can be fixed by enforcing that the quark susceptibility 
for a light quark with mass m in Eq. I|18l) : 



d(uu) 



dm 



(20) 



has a peak at the desired (/is,T) point. This has been done for the most accurately determined thermodynamical 
parameters, which are those in Pb-Pb collisions at v^attv = 

17.3 GeV, yielding T = 170 MeV (see Fig.©. Thereby, 
all relevant parameters are now fixed. The only free parameter left is the single cluster proper volume V c , which 
may vary as a function of center-of-mass energy and colliding system. The dependence of As on this parameter is a 
crucial issue to be studied because it can possibly establish a relationship between color deconfinement and strangeness 
enhancement. 



A. Heavy Ions 

In HIC, calculations are much simpler as the total volumes (namely the sum of proper volumes of all the clusters) 
turn out to be so large that the effect of the canonical suppression related to exact conservation of flavor quantum 
numbers can be disregarded 01 1 provided that flavor quantum numbers are distributed among clusters according to 
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the function JJjJ . Indeed, it should be pointed out that different scenarios have been devised in which the volume within 
which strangeness exactly vanishes (strangeness correlation volume) is less than the total volume, but this approach 
cannot account for the low yield of cj> meson [^j , so we will stick to the identification between the strangeness 
correlation volume and the total volume. If V c is large enough, quark multiplicities can thus be estimated by means 
of Eq. (|17|) . the grand-canonical formula, and As turns out to be independent of V. In fact, it must be emphasized 
that the total volumes fitted in the SHM for HIC, though affected by large errors, lie in a range (see Table where 
canonical suppression is negligible and As is sensitive to V c if this is roug hly below 10 fm 3 (see Fig. 0J. This little 
sensitivity to V c makes As certainly not a clear-cut probe for deconfinement, though hadronizing clusters as small as 
some fm 3 can indeed be excluded. 

Hence, by taking a full grand-canonical approach, the cross-over line for chiral symmetry breaking has been calcu- 
lated by minimizing the light quark mass susceptibility Xm- As has been mentioned, by forcing the location of the 
minimum to coincide with the fitted T and \xb in Pb-Pb collisions at y/s NN = 17.3 GeV, we have been able to set 
the parameter To to 170 MeV (see Fig- EJ- F° r this calculation, the u and d chemical potentials have been assumed 
to be equal, so that \u and Xd coincide with x-m, whereas the strange quark chemical potential /i s has been set to 
zero due to S = constraint. It must be noted that, in principle, \x u and fid differ because: 



, 2 

A*« = — + g^Q 
fi B 1 

Md = -jj T^Q 

= — - g^Q - fJ-s (21) 

where /xq is the electrical chemical potential, which vanishes only in isospin symmetric systems. However, this term 
turns out to be negligible for all examined collisions (see TableQJ and we have simply taken /i u = \i& = A*b/3. Having 
set To, the cross-over curve can be predicted from the model and this is shown in Fig. [5] along with fitted SHM points 
which are in satisfactory agreement with the calculation. 

It must be pointed out that in a grand-canonical finite system with fixed volume, the conservation of a definite 
initial values of baryon number, electric charge and strangeness necessarily leads to different values of the chemical 
potentials for the two different phases, such as the considered constituent quark phase and the hadron gas. Certainly, 
the volume is not supposed to be fixed in the transition, thus the assumption of equal /is's may be correct. On 
the other hand, once \xb has been set, the constraints S — and Q/B = Z/A lead to definite /j,q and us (beware 
the difference with ji s ) values which, in principle, are different in the two phases. While retaining /i s — 0, we have 
checked the accuracy of the assumption [aq = by studying the effect on As of non-vanishing electrical chemical 
potentials in the quark phase. The procedure is as follows: first, we have calculated As and quark masses with the 
main assumption, i.e. fj,Q = 0; then, by using those masses, we have calculated /xq in the quark phase by enforcing 
Q/B — Z/A; finally, with the obtained (1q, we have recalculated quark masses and the new As- The difference 
between the As's calculated in the two ways ranges from 1% in Pb-Pb collisions at SPS energy up to 2.7% in Au-Au 
collisions at AGS energy and is therefore negligible throughout. 

The parameter As has been calculated along the cross-over curve by using Eq. (|T7Jl and the comparison with the 
SHM fitted values is shown in Fig. [BJ The main prediction of the NJL-based model is an increase of As for decreasing 
yfs which is driven by the increase of /is (see Table HJl . This involves an enhancement of relative strange quark 
production with respect to u, d, the so-called Pauli blocking effect. On the other hand, the observed As does not keep 
growing but undergoes a dramatic drop at very small energies (see Fig. QJ, a fact which can be possibly explained 
by the onset of a purely hadronic production mechanism at low (yet not easy to locate) center-of-mass energies with 
possible local strangeness conservation |3^| . Looking at the deviations of the SHM-fitted central values of As in Au-Au 
and Si-Au from the theoretical curve in Fig. JSJl, it can be argued that this hadronic production mechanism takes 
over at energies as low as AGS's. Certainly, error bars are large and a definite conclusion needs more precise data. 
It should also be pointed out that a mild strange canonical suppression sets in already in AGS Si-Au collisions [To| 
which is not taken into account here. Another sizeable discrepancy between data and model shows up in Au-Au 
collisions at RHIC (y/s — 130 GeV). However, in this case, As has been determined by fitting ratios of hadronic yields 
measured at mid-rapidity, unlike all other quoted collisions for which full phase space multiplicities have been used. 
This method may have led to an overestimation of As if the midrapidity region is enriched in strangeness as observed 
in Pb-Pb collisions at SPS [H- 

It is now worth discussing the robustness of these results in some detail. The calculation of As based on Eq. H17|) 
depends on constituent quark masses and the UV cutoff parameter A besides the thermodynamical parameters T and 
fiB which have been set by using hadronic fits. This expression of As is indeed a general one for effective models with 
four fermion interactions in the mean-field approximation. It is then worth studying As by taking the constituent 
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masses as free parameters instead of fixing them by means of a particular effective model as we have done so far. 
Whilst M Ui d(fi, T) is strongly dependent on T , M s (fj,, T) is fairly independent of it (see Fig.|2J). Therefore, as far as the 
finite temperature sector is concerned, one can fairly conclude that this NJL model has actually one free parameter, 
either To or M Ut d(pBc,T c ). In other words, fixing To amounts to set a definite value of the light quark mass at the 
critical point. 

However, the resulting As value does not strongly depend on this parameter either, as shown in Fig. where the 
region in the (M U ,M S ) plane (with both A = 631.4 MeV and A = +oo) allowed by the fitted As values in Pb-Pb 
collisions at \fs NN = 17.3 GeV is shown: the bands pattern clearly indicates that As has a strong dependence on 
M s and much milder on M u ^. Furthermore, the M s value at the critical temperature is quite constrained by its 
corresponding value at zero temperature and does not undergo strong variations as a function of the temperature and 
To, as shown in Fig. (J2J, which is indeed a quite general feature of effective models |22| . Summarizing, we can state 
that this NJL model yields a As value in good agreement with the data (as determined through the SHM) essentially 
because it has a fit value of the constituent strange quark mass at zero temperature and density (about 500 MeV), 
whilst the particular value of To has much less impact on it (see also Fig. J2Jl). On the other hand, for fixed constituent 
quark masses, As is sizeably affected by the cutoff A. The calculation with A = +oo is meant to give some upper 
bound on its variation due to neglected UV contributions, though a re-analysis within a renormalizable model would 
be desirable to have a more accurate estimation. 



B. Elementary collisions 

In EC the fitted total volumes are small and the calculations have to be carried out within the canonical formalism. 
As has been mentioned in the previous section, the partition function to be used is that in Eq. (|12|) with constituent 
masses replacing current masses and a A cut-off as upper bound in the momentum integration in Eq. Unlike in 
the grand-canonical ensemble, As depends on the volume and a thorough comparison with SHM fitted parameters is 
much more involved. This is true in many respects: first, the fitted total volumes in the SHM are subject to large 
errors and obtained for point-like hadrons, thus reasonably underestimated. Secondly, for finite volumes one should 
in principle refit To so as to obtain a cross-over point for a definite total volume V and cluster volume V c , for each 
process, at the desired temperature. Finally, one should minimize the full expression of free energy F = —TlogZ 
(with Z as in $12^ ) to determine the constituent masses at finite volume. For the present, this calculation is not 
affordable as the minimization of functions involving five-dimensional numerical integration as in Eq. (|12[1 at each 
step with the due accuracy, implies exceedingly high computing times with presently available computers. For this 
reason, we have carried out a simpler calculation, with To kept to 170 MeV and constituent quark masses calculated 
in the model at the thermodynamic limit. 

Yet, even in this approximated calculation, the effect of flavor and color conservation over finite volumes can be 
studied and, hopefully, it can be verified if this mechanism of canonical suppression actually implies a reduction of As 
with respect to the thermodynamic limit consistently with the data. In fact, any conservation law enforced on finite 
volumes imply a reduction of heavy charged particle multiplicity stronger than lighter particles' so that a decrease of 
As for decreasing volumes is expected. The ultimate reason of this effect is the reduced energy expense needed, in a 
finite system, to compensate any charge unbalance with light particles in comparison with heavy particles. 

We can see this mechanism at play in Fig. [8] where As is plotted as a function of the total volume V (over which 
flavor is exactly conserved), for different (color-singlet) cluster volumes V c for an initially completely neutral system, 
such as e + c~ . The temperature has been set to T = 160 MeV, which is a fair average of the SHM fitted values 
@>0,E1 m e + e~ and the constituent quark masses have been set to M s — 452 MeV and M u ^ — 112 MeV, which 
are the thermodynamic limit values. Interestingly, As features a fair stability over a large range of total volumes 
greater than « 20 fm 3 for V c < 10 fm 3 within the data band, with a mild increase up to an asymptotic value which 
turns out to be closer to the thermodynamic limit (yet not equal to it, see also Fig. (0}) for larger V c . As the volumes 
fitted in a point-like hadron gas model in e + e~ and pp collisions turn out to be larger than 15 fm 3 , but not larger 
than ~ 80 fm 3 , we can argue that a suitable V c can account for the observed values and the stability of As over 
a reasonably large volume range, taking into account that the actual volumes are certainly larger than their point- 
like estimates. Moreover, the mechanism of local color neutrality induces the reduction of As value with respect to 
the thermodynamic limit of ~ 0.31 which is needed to quantitative reproduce the data. Within this approach, the 
difference between As in high energy HIC and EC would be mainly the result of the non vanishing initial baryon 
density (from ~ 0.31 to ~ 0.45) and of local color neutrality in EC over small regions (5 — 10 fm 3 ). 

The situation is somewhat different for pp collisions, due to non vanishing initial baryon number and electric 
charge. In this system, one can observe two compensating effects, namely Pauli exclusion principle favoring relative 
ss production and canonical suppression, favoring relative uu and dd production. As a result, As is nearly constant 
around 0.3, in apparent disagreement with the observed value which is about 0.2 (see Fig. EJ). The reduction of As 
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brought about by the decrease of V c is not sufficient to restore the agreement even for very low V c . The disagreement 
can be possibly explained by the inadequacy of the taken assumptions, such as the statistical distribution of charges 
among clusters according to Eq. JJJJ. 

V. CONCLUSIONS AND OUTLOOK 

We have presented a statistical model to explain the observed pattern of strangeness production in both elementary 
and heavy ion collisions. The basic idea is that full chemical equilibrium is locally achieved at the level of constituent 
quark degrees of freedom (within the framework of a simple effective model), at a different stage of the evolution 
process in elementary collisions (late) with respect to heavy ion collisions (early). In this approach, hadron formation 
takes place through the coalescence of constituent quarks and this accounts for the known observation of an incomplete 
strangeness equilibrium at hadron level in the statistical hadronization models. The underlying assumption is that at 
least the ratio s/u survives in the final hadrons. 

Besides the effect of different initial excess of matter over antimatter, the smaller relative strangeness production 
in EC with respect to HIC has been related to the smaller overall system size and to color confinement over much 
smaller regions. Particularly, the existence of a characteristic cluster volume, roughly between 5 and 10 fm 3 and 
independent of center-of-mass energy, can possibly account for the observed stability of As, along with the constancy 
of temperature and the weak dependence of As itself on the total volume V (if not below ^60 fm 3 ), at least in 
e + e~ and pp collisions. However, the dependence of As on cluster volume in HIC is so mild that it is not possible to 
use strangeness production to prove the formation of large color-neutral regions, i.e. color deconfincmcnt. 

The numerical analysis, carried out within an effective Nambu-Jona-Lasinio model, involved only one parameter to 
be adjusted which has been used to get the critical point agreeing with that fitted in Pb-Pb collisions. A satisfactory 
agreement with the data has been found in several heavy ion collisions and e + e~ , pp collisions. Deviations have been 
found at low energy HIC, such as Au-Au and Si-Au at AGS, where the strangeness production mechanism could be 
predominantly driven by hadronic inelastic collisions. Furthermore, a significant discrepancy has been found in pp, 
which might be perhaps cured by taking a different scheme of quantum numbers distribution among the hadronizing 
clusters. It should be pointed out that the analysis in EC has been performed with some approximations as a full 
consistent treatment of the Nambu-Jona-Lasinio model for finite volumes is presently beyond our possibilities. 

The model-dependence of our results has been discussed. We are reasonably confident that the obtained results are 
quite general and stable within non-renormalizable models with four-fermion interactions treated in the mean field 
approximation. The exploration of renormalizable effective model is under investigation. 

APPENDIX A 

We want to prove that: 

N 

E Il^Q^Q^Q. ( 22 ) 
Qi,...,Qjv i=l 

is the global partition function in Eq. I|12|) . with Zi as in Eq. 0, provided that all clusters have the same volume V c . 
We will first prove it in the canonical framework, with each cluster having also the same temperature T . One can 
rewrite Eq. I|22[) by using Eq. (JTJ and the integral representation of the Sq^^ as: 

e iQi-<& J ' dn{6 u Q 2 ) eWMa,^) (23) 

where Z^ has been written in the simple form exp[V c f(9i, 6*2, 4>i)} according to Eq. ©. It must be stressed that this 
expression of Z is very general and does not depend on the Hamiltonian of the considered model. Moreover, the 
function / is the same for all clusters for it depends only on the temperature which has been assumed to be the same. 
We can now perform the summation over all Qi in Eq. ill' 21) and get: 

N 

X- giSiQi • (0, -4>) = J]( 27 r) 3 5 3 {4> - </>,) (24) 

Qi,...,Qjv i=l 



E 



n 



2tt 



= iQ ■ - lEiQi • <fi 



N 

n 



n 



~2T 
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so that the integration over <j)ji in Eq. (|22() can be easily performed and one is left with: 



n 

i=i 



2tt 



n 



y c f(0i,e 2 ,<i>) 



The integral over SU(3) group is the same for all clusters and Eq. (|25|) can then be written as: 



3 



(1^,02) z i {e 1 ,e 2 ,<j>) 



where V = T^Vi = NV C , which is equal to the right hand side of Eq. 1)12(1. 



N 



v/v c 



(25) 



(26) 



APPENDIX B 



We show that the configuration probabilities u>(Qi, . . . , Qjv) in Eq- @ minimize the free energy of a system with 
volume V and temperature T which is split into N color-singlet clusters with volumes V\, . . . , Vat such that ^ = V . 
Let p be the full probability of a single state in this system and w the probability of a configuration (Qi, . . . , Qjv). 
Then: 



jn n \ TT exp(-£ , i/T) 

Pstate = . . . , t^Arj — ^ OQ ; ,Q i;3 Osinglcti 



(27) 



Li 

where Zi is given by Eq. J7J), Ei is the energy of a cluster, Qi :s is the vector of flavor quantum numbers of the state 
and ^ s i n gicti signifies the color singlet constraint on each cluster. The entropy reads: 

s = - E p 1o &p 



a —Ei/T 



Ql,...,QiV statesi statcsjsj 



logw(Qi,...,Qjv) 



JV 

n 



exp(-Ej/T) 



(28) 



where the flavor and color constraint on each cluster arc now implied in the sum over the states. 
Eq. (|2~5j) can be worked out: 



S 



n —E-/T 



Qi,...,Qiv states i states^ 

E, 



N 



log w(Qi, . . • , QaO - - log l[ Zi(Qi) 



i=i 



= t«(Qi,..-,Q 

Qi,..,Q]v 



JVj 



logw(Qi 



x, ■ ■ ■ , Qjv) + ^ l(gi) j qi W + logfl Zi(Q 



(29) 

where 2i(-^»)(Qi,...,Qjv) ^ s the average energy for the i th cluster and a fixed configuration (Qi, . . . , Qjv)- The total 
average energy of the system, to be identified with the internal energy U, evidently reads: 



jY 



(Qi,...,Qjv) 

Qi,— ,Qn *=l 



(30) 
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so that entropy can be rewritten as: 



s = E -f(Qi,.»,Qw)logw(Qi,...,Q 

Qi,...,Qjv 



T 



^ «;(Qi,...,Q w )logJJZ i (Qi) 



(31) 



Qi,...,Q, 



and the free energy F as: 



F = U-TS = T J2 w(Qi,...,Q w ) 
Qi,...,Qiv 



JY 



Iogi£;(Q 1 ,...,Q^)-bgJJZ i (Q i ) 



(32) 



Now we seek for the probabilities w which minimize the free energy with the constraint Y]w = 1. Therefore we have 
to introduce a Lagrange multiplier A and look for the extremals of F + \(^2 w— 1) with respect to each u>(Qi, . . . , Qw) 
and A. This leads to: 



/v 



log«;(Qi,...,Qj V ) = logJJZ i (Q i ) 



T 



(33) 



and 



after proper normalization. 



w(Qi, . . . , Qjy) - — i-rjv ~ , u 

L,Q U ...,Q N lli=l A(,*4i)0Q,EiQ» 
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TABLE I: Parameters determined through statistical model fits to measured multiplicities in some elementary and heavy ion 
collisions. Heavy ion parameters have been taken from refs. |sl l3^l whilst elementary collisions from refs. p3. ll2f : e + e~ at \fs = 
91.2 and pp points have been refitted with updated data and hadronic input parameters. The pp point has been fitted with a 
different parametrization of the strangeness suppression |T^ | and 7s has been replaced with the mean value of strange quark 
pairs (ss); the corresponding 7s would be ~ 0.5. The numerical values in Pb-Pb at y/H = 8.7 GeV have been obtained in 
ref. |36| with preliminary data from exp eriment NA49, those in Au-Au collisions at y/s = 130 GeV from a numerical analysis of 
midrapidity ratios as quoted in ref. |37|l . The quoted central values of volumes (the volume is defined as the sum of the volumes 
of all produced clusters), have been obtained with pointlike hadrons, thus they must be taken as lower limits for the actual 
values. The errors on volumes, not quoted here, are generally large and can be up to 50%. Also quoted the central values of 
the electrical chemical potentials obtained in heavy ion collisions. 



Collision - 


yi (GeV) 


T (MeV) 


fj-B (MeV) 


7s 


As 


HQ (MeV) V (fm 3 ) 


Au-Au 


130 


167±7.2 


45.8±6.4 


1.04±0.10 


0.476±0.049 


-1.42 




Pb-Pb 


17.3 


158.1±3.2 


238±13 


0.789±0.0582 


0.447±0.025 


-6.87 


3460 


Pb-Pb 


8.7 


149.0±2.4 


393.7±8.3 


0.822±0.058 


0.585±0.052 


-11.2 


2067 


Si-Au 


5.4 


133.4±4.3 


581±32 


0.845±0.101 


0.72±0.14 


-10.7 


330 


Au-Au 


4.8 


121.2±4.9 


559±16 


0.697±0.091 


0.43±0.10 


-12.4 


2805 


e + e~ 


14 


167.4±6.5 




0.795±0.088 


0.243±0.036 




15.9 


e + e~ 


91 


159.2±0.8 




0.664±0.014 


0.225±0.004 




52.4 


PP 


27.4 


162.4±1.6 




(ss) 0.653±0.017 0.201±0.005 




25.5 


PP 


200 


175±11 




0.491±0.056 


0.214±0.025 




35.5 


PP 


900 


167.±9.0 




0.533±0.054 


0.230±0.033 




77.3 



FIG. 1: As in various elementary and heavy ion collisions (from 36]). 
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FIG. 2: Constituent masses of the quark s and the light quarks u, d as a function of temperature, with fixed fJb u ,d/T = 0.5 
and jj, s — 0. The three curves correspond to different To values. The vertical lines define the la band for Pb-Pb collisions at 



fs NN = 17.3 GeV, as fitted in the SHM 9] 
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FIG. 3: Chiral susceptibility Xm for the quark u or d with current mass 5.5 MeV as a function of temperature T, for different 
values of the To parameter in the NJL model. The baryon chemical potential fiB has been set to 238 MeV and fiQ has been 
set to 0, in accordance with the analysis of full phase space Pb-Pb data at SPS energy 0. The maximum of — Xm is located 
at T = 158 MeV, the fitted value in Pb-Pb collisions at SPS energy for T = 170 MeV. 
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FIG. 4: Calculated As in the NJL model in a color singlet cluster with T = 158 MeV, baryon chemical potential /kb/T = 1.5, 
Hq = and fi s = 0, as a function of the volume. The dashed line indicates the grand-canonical value of As. The constituent 
quark masses have been fixed to their values in the grand-canonical limit: M u = Md = 103 MeV and M s = 452 MeV. 
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FIG. 5: Cross-over curve in the NJL model with To — 170 MeV in the \{ub,T) plane. Black dots with error bars are the values 
obtained within the SHM analysis in various heavy ion reactions 
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FIG. 6: As calculated along the cross-over curve predicted in the NJL model with To = 170 MeV. Black dots with error bars 
are the values obtained within the SHM analysis in various heavy ion reactions 0, [3(J . 
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FIG. 7: Allowed regions of constituent quark masses M Ul di M 3 determined by the parameters As = 0.447±0.025, T = 158.1±2.3 
and fiB = 238 ± 13 MeV as fitted in Pb-Pb collisions at i/s NN = 17.3 GeV with /j,q — and /i s = 0. The lighter hatched 
region corresponds to the A momentum cut-off as in the NJL model [l9|l whereas the hatched darker region to the free constituent 
quark gas (i.e. A = oo). The solid line shows the predictions of the full NJL model with To ranging from zero (no KMT term) 
to oo (CASE I in ref. [13 ]') from left to right, in the grand-canonical ensemble with T = 158.1 MeV and fiB = 238 MeV while 
the black dot shows the values of masses for To = 170 MeV. The dash-dotted horizontal and vertical lines correspond to the 
T = fj,B — constituent masses values for the strange quark and for the light quarks u,d respectively. 
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FIG. 8: Plot of As as a function of the total volume V of the hadronizing clusters, for various (color singlet) single cluster 
volumes V c in an initially completely neutral system. The constituent quark masses have been set to their values obtained 
in the grand-canonical case with T = 160 MeV and all chemical potentials vanishing. The dotted line has been obtained by 
disregarding the color singlet constraint and keeping only the flavor constraint, for a given volume. The arrow indicates the 
asymptotic value in the grand-canonical limit. The horizontal band is the region determined by the maximal spread of central 
values obtained in the SHM fits in e + e~ |3.ll2LlT^| and in pp collisions Q|. 
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FIG. 9: Plot of As as a function of the total volume V of the hadronizing clusters, for various (color singlet) single cluster 
volumes V c in an initially pp-like system (B = Q — 2, S — 0). The constituent quark masses have been set to their values 
obtained in the grand-canonical case with T = 160 MeV and all chemical potentials vanishing. The dotted line has been 
obtained by disregarding the color singlet constraint and keeping only the flavor constraint, for a given volume. The arrow 
indicates the asymptotic value in the grand-canonical limit. The horizontal band is the region determined by the maximal 
spread of central values obtained in the SHM fits in pp collisions [2lll2|. 
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